Thermal Mass Spectra of Scalar and Pseudo- Scalar Mesons in 
IR-improved Soft- Wall AdS/QCD Model with Finite Chemical 

Potential 

Ling-Xiao Cufl and Yue-Liang WiJll 

State Key Laboratory of Theoretical Physics(SKLTP) 
Kavli Institute for Theoretical Physics China (KITPC) 

Institute of Theoretical Physics, UCAS 
Chinese Academy of Sciences, Beijing, 100190, China 
(Dated: February 21, 2013) 

Abstract 

The thermal mass spectra of scalar and pseudo scalar mesons are studied in an IR-improved soft- 
wall AdS/QCD model. The Reissner-Nordstrom AdS black hole metric is introduced to describe 
both the finite temperature and density effects. The thermal spectral function is computed, which 
shows that the position of peak moves towards a smaller value and the width of peak also broadens 
to a larger value as the temperature increases. The critical temperature at which the peaks are 
completely dissolved has been found to be around Tc — 160 ~ 200 MeV. The effect of chemical 
potential is shown to be the same as the one caused by the temperature. It is found that when the 
temperature approaches to zero T ~ 0, the melting down of mesons occurs at the scaled critical 
chemical potential /Uq/k ~ 2.2 ~ 2.4 GeV. 

PACS numbers: 



*Elcctronic address: 'clxyx @itp.ac.cn| 
''"Electronic address: |ylwu@itp. ac.cn 



1 



I. INTRODUCTION 



The properties of Thermodynamics of quantum chromodynamics (QCD), which is related 
to the experiment at the Relativistic Heavy Ion Colhder (RHIC) and the Large Hadron Col- 
hder (LHC), has been attracted a great deal of attention Recently. The strongly correlated 
quark-gluon plasma (QGP) produced from RHIC involves many non-perturbative problems 
of QCD, such as chiral symmetry breaking and restoration, deconfinement phase transition, 
QCD phase diagram and so on. Many theoretical tools have been developed to describe 
those phenomena: effective QCD theories such as effective chiral dynamical model (CDM) 
Q% of low energy QCD, QCD Sum Rules and lattice QCD. Although the effective CDM 
can provide well prediction for the mass spectra of ground state mesons, while it is hard to 
characterize the linear trajectory behavior for the resonance meson states. When it comes 
to the dense matter situation, the lattice QCD calculation encounters difficulty for the sign 
problem in nonzero chemical potential case. 

In the framework of the gauge/gravity duality approach 046|, the holographic QCD has 
been developed to bring new light on solving the problem of nonperturbative QCD. Searching 
the bulk gravity theory which corresponds to realistic QCD was first pursued in [sf. Mainly 
there are two complementary ways to follow: the top-down approach and the bottom-up 
approach. The former starts with the string solutions and find the brane configuration for the 



dual gravity to reproduce some basic QCD features, e.g. D3/D7 [7H9j, DA/D8/D8 |lOl 111 
D4/ DQ/ D6 [12}. The latter known as a phenomenological AdS/QCD model consists of a 
gauge theory in a curved space with some field contents which are dual to some operators 
in QCD. The so-called hard- wall AdS/QCD model with a sharp infrared (IR) cut off AdS^ 



metric was developed firstly in Refs. |13Nl6|. It can correctly realize the chiral symmetry 



breaking and low-lying mesons states, while the higher excited states of hadrons were found 



to deviate the Regge trajectory. Then the soft-wall AdS/QCD model [17| was constructed 
to remedy this drawback by introducing a background dilaton field which suppresses the 
action gradually at IR region. Applying the WKB approximation analysis, it was shown 
that the dilaton's quadratic behavior at IR region will lead to the Regge trajectory for the 
excited meson spectra. However, the chiral condensation in the soft-wall model was found 
to be proportional to quark mass. A large amount of works have been done to incorporate 
both the chiral symmetry breaking and the linear trajectory of excited meson spectra in the 
AdS/QCD model. In [l8|, a quartic interaction term has been introduced into the action 
in soft-wall model to realize the explicit and spontaneous chiral symmetry breaking, while 
such a term was shown to cause an instability of the scalar potential and a negative mass for 



the lowest lying scalar meson state. In Refs. |19l . |20| , by simply modifying the bulk gravity 
at IR region, a consistent soft- wall AdS/QCD model has been constructed and the resulting 
mass spectra for light resonance states have been found to agree with the experimental data. 

It is also interesting to investigate the finite temperature and density properties in the 
Ad S/Q CD model. The hard- wall's finite temperature effects was studied in j2l[. In Ref. 

;he scalar glueball and light scalar mesons spectral have been studied in the soft- 
wall model and the highest temperature for which the meson state can exist was found to be 
about 40 — 60 MeV. Compared with the Hawking-Page (HP) transition temperature in the 
soft-wall model, which was found to be Thp — 192 MeV in |25|, the meson states dissociation 
occurs in the confined QCD phase, which is far from the deconfinement transition. One way 
to increase the temperature is to set up a different mass scale in the soft-wall model. In 



Ref. |26|, |27|, J/ip thermal spectral function was studied by treating the scale as a flavor- 
dependent parameter. In Ref. j28|, the melting of J/ip spectral peak was found to occur 
at T !=a 540 MeV. However, the light mesons' dissociation temperature is still too low. 



The chemical potential is another interesting properties, which has been widely investigated 



in holographic QCD. In Ref. |29H31|. the spectral functions have been studied by taking 
D3/D7 setup with fundamental matters at finite baryon density. Density effects on the 
spectral function in the soft-wall model have been discussed in |32l . 

In js^l , we extended the the IR-improved soft-wall [l9l to a thermodynamic model at fi- 
nite temperature. Spectral functions of vector and axial-vector mesons have been computed. 
The highest temperature for which the meson state can survive was found to be around 
Tc ~ 200 MeV without quartic interaction term. However, difficulty remains in obtaining 
the spectral function of pseudo-scalar meson due to the coupling between the axial-vector 
field and the pseudo-scalar field. In this paper, we only pay attention to the pseudoscalar and 
scalar mesons. To be completeness, we will construct the soft- wall AdS/QCD model which 
works well in both zero temperature and finite temperature region for the pseudoscalar and 
scalar mesons. In Sec. [Tll inspired by Refs. 



19, 20 



we will build a thermodynamic model 
for scalar part by introducing a modified Reissner-Nordstrom AdS black hole. In Sec. IIIH 
parameters are fixed at zero temperature with density. The effects of finite temperature and 
density will be discussed in Sec. IIVI Our conclusions and remarks are presented in the final 
section. 



II. IMPROVED SOFT-WALL ADS/QCD MODEL 



In this section, we will construct a phenomenological thermodynamic AdS/QCD model 



inspired by the predictive holographic QCD model [19|, |20[. In Ref. |34| . we studied the 
spectral functions of the vector and axial- vector mesons. However, the pseudo-scalar part is 
the most intriguing one due to the coupling with axial-vector field. As our main goal is to 
investigate the finite temperature effects, we will consider the simple case by taking off the 
gauge field part and focus only on scalar and pseudo-scalar part in the 5D action. Thus the 
5D action becomes: 



S = j d^x V^e-*(^) Tr [\DmX\ 



m 



X 



\X\ 



(1) 



with m'j^ 



3. The quartic term A|X|^ is introduced to improve the ground state of scalar 



mesons masses which have been discussed in Ref. |19|. $ is the dilaton background field. 
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which can lead to a linear trajectory if the UV behavior is taken as ^{z — )■ oo) ~ z' 
The fundamental field X, which corresponds to a dimension three scalar operator q%q^, can 
be written as: X{x,z) = {X(z) + S{x, z)e^^'^^^'^^) . S{x,z) and 7r{x,z) are scalar field and 
pseudo-scalar field respectively, X{z) = {X{x,z)) is the bulk vacuum expectation value. 

In order to investigate finite temperature and density effects, we need to find a dual geome- 
try in the bulk. According to AdS/CFT correspondence, the quark chemical potential, which 
is the coefficient of the quark number operator qj^q, corresponds to the time-component of 
U{1) gauge field in the bulk. The gravity action which describes the interaction of the U{1) 
gauge field with the AdS space is given by: 



S 



1 



(2) 



When using the following ansatz: 



At = Atiz), 

A = o 



1, 



dz' 



which leads to a well-known solution called Reissner-Nordstrom AdS black hole: 

4 / \ 6 



AAz) 



Q 2 



(3) 
(4) 

(5) 



(6) 
(7) 



with ^ Q ^ a/2, which is the charge of the gauge field. 

In this paper, we will adopt a modified Reissner-Nordstrom AdS black hole inspired by 
the one introduced in [3^ 



ds^ = a\z) {f{z)de-^dx\ 



i=l 



dz^ 



(8) 



In order to incorporate both chiral symmetry breaking and linear confinement, the IR region 
of the metric is modified as: 

a\z) = l/z' + fil (9) 

Here fig is a constant mass scale which characterizes QCD confinement. As shown in (lol. |20|. 
such an improved 5D metric will lead to a non-trivial dilation solution. f{z) is the same 
as the one in Eq. [TJ so the Hawking Temperature, which corresponds to the temperature in 
boundary theory, is not changed under this modification: 



H 



1 


df 






An 


dz 







(10) 



The chemical potential can be defined by the condition that gauge field vanishes near the 
black hole horizon: At (zh) = 0, so we have 



Q 

jj, = K 

Zh 



(11) 



At is a free parameter, which can be determined in many different ways [35|, |36[. In this 
paper, we will directly give the result of h/k. 

The bulk vacuum expectation value of the fundamental field X has the following form 
for the two flavor case 



{X{Z)) = -V{Z) l2. 



(12) 



While the VEV v{z) is related to the dilaton field ^{z) through the equation of motion: 



Q = dAa^[z)f{z)e-''^^^dM'z] 



m 



X 



- I V[Z] 



(13) 



In this paper, we take the same form of VEV v{z) proposed in [3J]. It consists of two 
parts: zero temperature part Vq{z) and the finite temperature part Vi{z) In f{z). 



v{z) = vi{z)lnf{z) +vo{z), 



(14) 



where the zero temperature part vq{z) is taken from the type lib model in Ref.[l9[ which is 
given by Table. HI Its IR boundary behavior vq{z — )■ oo) ~ -y/z can lead to a more accurate 
prediction for all the mass spectra in zero temperature and density. 



Model 



^'o(^) 



Parameters 



lib z{A + Bz'^){l + Cz^)-^/^ A 



m,C, B = a/C, C = (Byifi,^^))^/^ 



TABLE I: lib model which is taken from [3] with = [idl\f^. The numerical values are shown 
in Tabel. IE 



The finite temperature part is introduced to regulate the singular behavior of $(-z) near 
the horizon 3J] . From Eq. [13] , we can obtain the non-trival form of the dilaton. In the 
vicinity of horizon, dilaton field $(2;) has a divergent behavior: 



ln/(^) 

vx{z)f'{z^ 



-T{z) + 



(15) 



with 



T{z)^?>a\z)v^{z)\j\z)v^{z) (16) 

Note that near the horizon, fi^z^) — )■ 0, the term In /(z) becomes divergent. If v\{^z^ takes the 
form as the solution of Ti^z) = 0, the dilaton gets smooth solution at the horizon as shown 
in js^]. In this work, Vi{z) becomes too complicated due to the introduction of chemical 
potential. In zero density limit, ^1(2;) has the same form as in 341 . 



vi{z) = c„ exp 



nix f ^ 22 



(17) 



When considering the effects of the quartic interaction term, the coefficient of the quartic 
term A takes the following form 



A 



Ao 



cx{\nf{z)y 



with 



P>3, 



(18) 



For simplicity, we still take p = 3 in current consideration. Thus all parameters in this 
model will be inputted by experimental data at zero temperature and density limit in next 
section. 



III. MODEL AT ZERO TEMPERATURE AND DENSITY LIMIT 



A. Fitting Parameters 



We will make numerical calculations for the mass spectra of scalar and pseudo-scalar 
mesons in zero temperature and density limit. Turning off chemical potential and temper- 
ature, which means f{z) — 1 and Q ^ 0, the metric takes the following form which is the 



same as the one used in [19l 20 



ds^ = a'iz) [df-Y, dxl -dzA, (19) 



i=l 



The finite temperature part of VEV vanishes, left only the zero temperature part. 

v{z)=v^{z), (20) 

/ig as a mass scale can be fixed from a global fitting. The left three parameters m^, a and 
7 are fixed by the well measured experimental value of ttlt^ = 139.6 MeV, the Gell-Mann- 
Oakes-Renner relation /^m^ = 2mqa, and the mass of scalar mesons. The values of three 
fitting parameters are presented in Table. HTl 



Aq mq (MeV) a 5 (MeV) (MeV) 7 
Ao = 430 268 445 oT~ 

Ao = 15 6.99 228 505 0.124 



TABLE II: The numerical values of parameters for the case without the quartic term (Aq = 0) and 
the ones with the quartic term (Aq = 15), with C, = \/3/(27r) and fig = /i^/i/S- 



B. Mass Spectra of Scalar and Pseudo-Scalar Mesons 

The fundamental scalar field can be decomposed as X(x, z) = (v(z)/2 + S{x, z))e^^'^^^'^\ 
where S{x, z) is the scalar meson field and 7r(x, z) = 7r"(x, z)t°- the pseudo-scalar field. Now 
we can derive the equations of motion for both two fields: 

S : S:{z) + S:{z) - $'(z)) + S^{z) (ml + 3a(z)^ + ^Aa(z)2^(z)^) = 0(21) 

PS : n:iz) + ttK.) + ^ - ^'(^)) + ^l^niz) = 0, (22) 

(23) 

Under the boundary condition Sn{z — 0) = 0,dzSn{z 00) = 0;7in{z — )• 0) = 
^idzTTniz — 00) = 0, we numerically solve the equations by shooting method. In this 
way, we find out the normalizable modes for scalar and pseudo-scalar fields in the bulk and 
their eigenvalues, which are the masses of mesons. The predictive mass spectra for the scalar 
and pseudo-scalar mesons are showed in Table. IIIII Note that the scalar states /o(980 ± 10), 
/o(1505 ± 6), /o(2103 ± 8), /o(2314 ± 25) should be classified into the isosinglet resonance 
scalar states of SU (3) octet mesons, rather than the SU (3) singlet resonance scalar states. 
Here to circumvent the difficulty of the finite behavior of pseudoscalar mesons, we have ig- 
nored the coupling between pseudoscalar field and axial- vector field from the action [TJ This 
will cause the tiny mass splittings between scalar and pseudo-scalar meson mass eigenvalues, 
especially the ground state mesons. After adding the quartic interaction term, mass spectra 
for scalar meson can be improved a little. 
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experimental. (MeV) A = A = 15 


/o experimental. (MeV) A = A = 15 





139.6 


139.6 139.6 


"J 'J "-'-150 


139.8 275.3 
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1350 ± 100 


1439 1431 


1350 ± 150 


1439 1453 
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1816 ± 14 


1701 1738 


1724 ± 7 


1701 1758 


3 




1926 1967 


1992 ± 16 


1926 1987 


4 




2126 2159 


2189 ± 13 


2126 2179 


5 




2308 2326 




2476 2346 


5 




2475 2476 




2633 2495 



TABLE III: The experimental and the predicted mass spectra for Scalar (right side) and Pseudo- 
Scalar (left side) meson 

IV. MODEL WITH FINITE TEMPERATURE AND DENSITY 

When temperature and chemical potential are turned on, the equation of motion is given 
as follows in momentum space by performing the Fourier transformation(we put three- 
momentum to zero = 0): 



V a{z) f[z) J 



(!M£g£)! H- M£mfl±^) ^ 0. (24) 

PS : A-,., ^ (M£) + I M£) _ . ^ ^ 0. (26) 

V o.{z) f[z) v[z) J f[zY 

As the temperature and density increase, the horizon of black hole Zh moves from infinity 
to boundary side. Solutions of equations of motion drop into black hole before they vanish. 
Therefore we cannot use the method of finding eigenmode. Alternatively, we will investigate 
spectral function which is the imaginary part of the retarded Green's function. 

First we will check the boundary behavior of the solution. Near the origin, we can 
extract the asymptotic solution of Eq. |2l]and Eq. [251 For convenience, we change the radial 
coordinate z to dimensionless variable m as u = z/ Zh- The two linear independent solutions 
are given as follows: 



S : i^i ^ mVi , u^Yx 



PS : i^,^.J,f^y K.-.uyA'^] (26) 



\Zh J \ Zh 

Here, Ji and Yi are the first-kind Bessel function and second-kind Bessel function respec- 
tively. On the other hand, near the horizon, we take the in-falling boundary condition which 



corresponds to retarded Green's function [37 



K^^{1- u) ^^^) (27) 



The solutions of equation of motions can be expressed by the combination of the two inde- 
pendent asymptotic solutions. The coefficients A{u, q) and B{uj, q) will be fixed by the IR 
in-falling boundary condition: 

K{u) = A{u, q)Ki{co, q, u) + B{u, q)K2{u, q, u) ^ {1 - u)'^^) (28) 
The on shell action of scalar field part reduces to surface terms: 

e-''^'^f{z)a{z)-^S{p, z)d.S{p, z) , (29) 



Following the prescription in |37|, after substitute Eq. [26] into surface terms, we find that the 
spectral function which relates to the imaginary part of two point retarded Green function 
is proportional to the imaginary part of B{uj, q)/A{(jj, q). The results for pseudo-scalar part 
is the same. 

p{u, q) = --Im G{u, q) 9{co' - q') a Im^^ (30) 



Let us now make a numerical calculation. As we have discussed in [3^, c„ and c\ play the 
role as regulators to regulate the singular behavior near the horizon, which have no effect 
on finite temperature and density behavior. In principle, we should assign their values as 
small as possible, while it turns out to be more difficult if we adopt extremely small values 
in practical calculation. In Ref. jlsj, we have computed the quark number susceptibility at 
finite temperature by taking the limit that momentum and frequency are eventually taken 
to be zero. In this case, the modified terms which contain c„ and cx can be neglected. In 
this paper, we take = cx = 10"^. 



A. Critical Behavior with Finite Temperature { T ^ and i_i = 0) 

In Fig. [H we turn off chemical potential and check the finite temperature behavior in 
scalar and pseudo-scalar channel. As the temperature increases, the horizon Zh approaches 
boundary so that both scalar and pseudo-scalar mesons become unstable and melt. For the 
low temperature, the positions of the peaks are in accord with the masses which have been 
given in zero temperature case in Table. IIIII Increasing temperature, the peaks are shifted 
towards smaller values and the widths become broader. We can define a melting point, where 
no peak can be distinguished above this temperature. The temperature Tc ^ 200 MeV(A = 
0) or Tc ^ 160 MeV(A = 15) can be taken as the critical temperature of dissolution of 
mesons state. Here the critical temperature Tc ~ 200 MeV without the quartic term is 
consistent with other NJL models' predictions js], [39-42|, and also with the one obtained 



from the melting point for thermal mass spectra of the vector and axial- vector mesons [34 
. The value ^ 160 MeV with the quartic term approaches to the critical temperature 
yielded from the quark number susceptibility when the phase changes from the hadron phase 
to the quark-gluon-plasma phase at zero chemical potential j38[, and also to the one from 
the lattice QCD calculations j43|. 



We can get more information by fitting the spectral function with a Breit-Wigner form: 



(31) 




FIG. 1: The results of spectral function for scalar meson (left side) and pseudo-scalar meson (right 
side) . 




T (MeV) 



FIG. 2: The relation between the location of the first peak and the temperature for both scalar 
and pseudo-scalar parts with quartic interaction term. 

where m and F correspond to the location and width of the peak. P{uj'^) is representing a 
continuum which is taken the form -P(w^) = Ci + C20J^ + c^^u^Y*. To see the shift of the 
mass of resonance in both scalar and pseudo-scalar channels, we consider the lowest lying 
state with quartic term and plot the relation between the mass m and the temperature T 
in Fig. H 

It can be seen from Fig. [2] that in the low temperature region the location of the peak 
drops linearly as temperature increases. When temperature is near T = 100 ~ 150 MeV, 
the mass starts to decrease drastically, and at T = 150 MeV the mass is reduced to about 
30% of its value at T = 0. The lower the temperature is, the more difficult the numerical 
analysis is. This is because the horizon Zh becomes too large and the method cannot be 
used anymore at very low temperature. We have to adopt the way used in zero temperature 
calculation to find the eigenvalue of the solution. In this work, we shall not investigate such 
a region intensively. Just from the extrapolation of the fitted curve, we can find that around 



the critical temperature Tc — 160 MeV the mass seems to approach zero. The same results 
can also be found in some other model's prediction jsf, which has been shown to relate to 
the chiral symmetry restoration at finite temperature. 



B. Effect from Chemical Potential (T / and /U / 0) 

Let us now turn on chemical potential. From Eq. [TTl it can be seen that the effect of 
chemical potential is the same as the temperature and makes horizon Zh move close to the 
boundary side. In Fig. [3l we present the numerical results for the first peak of the scalar 
and pseudo-scalar mesons spectral function in different chemical potential. Here we have 
used the ratio /i/k instead of fitting the dimensionless parameter k. It is shown that the 
increasing of chemical potential leads to the same behavior as the temperature does. As 
the horizon moves towards boundary, the solutions are absorbed by the black hole and the 
bound state starts to collapse. 
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FIG. 3: The chemical potential effects on spectral functions for both scalar and pseudo-scalar 
mesons. 
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FIG. 4: Critical points of scalar and pseudo-scalar parts in T — fi planes 

As the increasing of both temperature and density can lead to an unstable state as 
expected, we can define a critical point where the first peak of ground state completely melt 
down and cannot be distinguished anymore. Finding out the critical values (Tc, /ic) for both 
temperature and density, we can plot it in the T — fi plane. The critical points for both 
scalar and pseudo-scalar parts are showed in Fig. HI When the temperature approaches to 
zero (T — )■ 0), the critical chemical potential is found to be Hq/k, = 2.4 GeV for scalar 



part and ^0/ n = 2.2 GeV for pseudo-scalar part. On the other hand, when the chemical 
potential goes to zero /i = 0, one obtains the critical temperature Tq ~ 200 MeV which is 
the same as the one given in previous subsection. The curve in the T — fi plane divides the 
whole plane into two regions. The region under the curve represents the phase in which the 
meson states still exist, which can be interpreted as the "hadronic phase". While the region 
above the curve corresponds to the real QCD's quark-gluon phase in which the meson states 
completely dissolve. As mesons are the bound states of quarks and antiquarks, the critical 
behavior of their disappearance can be understood as the chiral symmetry restoration. 

For the free dimensionless parameter k, one may get its value by comparing with other 
models' results. In jsof, there appears a relation To/fiQ ^ 1/2 for the critical chemical 
potential /iq at T = and critical temperature Tq at /i = 0, which then allows one to fix the 
parameter k ^ 1/6. 



V. CONCLUSIONS AND REMARKS 

We have investigated the IR-improved soft-wall AdS/QCD model with zero tempera- 
ture and finite temperature. An IR-improved Reissner-Nordstrom AdS black hole metric 



has been introduced to describe the finite temperature and density effects. Following |34 
an additional finite temperature part of bulk vacuum expectation value is introduced to 
get a smooth dilaton solution. In zero temperature case, the predictive resonance meson 
states agree well with the experimental data except for the ground state of scalar mesons 
which will further be investigated elsewhere. At finite temperature, we have computed the 
spectral function of mesons, a broadening of the peaks and moving towards the smaller 
values of the mass as the temperature increases have been demonstrated in detail. The 
critical temperature at which the peaks completely dissolved has been found to be around 
Tc ~ 200 MeV(without quartic term) and Tc ~ 160 MeV(with quartic term). By fitting the 
spectral function in terms of the Breit-Wigner form, we have quantitively studied the mass 
shift as the temperature increases. The feature that pion mass approaching to zero around 
the critical temperature agrees well with many other theoretical predictions. We have also 
investigated the effect of chemical potential. Its effect has been shown to be the same as the 
one caused by the temperature, namely it can lead to an unstable meson sate. It has been 
found that when the temperature goes to zero T ^ 0, the melting down of mesons occurs at 
the scaled chemical potential fio/n ~ 2.2 GeV for pseudo-scalar mesons and fj^o/n ~ 2.4 GeV 
for scalar mesons. 
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